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ABSTRACT 

From observations of a stream of gas bubbles 
rising through a liquid, a two-phase mathematical 
aodel Is proposed for calculating the induced turbu- 
lent vertical liquid flow. The bubbles provide a 
large buoyancy force and the associated drag an the 
liquid novas the liquid upward. The liquid pwped 
upward consists of the bubble wakes and the liquid 
brought into the Jet region by turbulent entrainment. 

The expansion of the gas bubbles as they rise through 
the liquid Is taken Into account. The continuity and 
MOMWntuM equations are solved numerically for an ax- 
laymetrlc air Jet submerged in water. Hater pumping 
rates are obtained as a function of sir flow rate and 
depth of submergence. Comparisons are made with lim- 
ited experimental information in the literature. 

NCNEH&ATUHE 

a local outer radius of liquid Jet region 

a c local radius of bubbly core region 
K q turbulent Jet entrainment coefficient 

g acceleration of gravity 

K ratio of bubble wake volume to bubble volume 

L depth of Jet origin below liquid surface 

M mass flow rate 

p pressure 

p ft atmospheric pressure 

Q volume flow rate 

r radial coordinate, origin is at Jet axis 

R perfect gas constant 

T absolute temperature 

u velocity in x direction 

u. terminal velocity of single bubble rising In 

Infinite liquid region « 

v velocity in radial direction 

Stoner faculty fellow, NASA Lewis Research Center. 


x vertical coordinate along Jet axis, origin is 
at gas release orifice 

p fluid density 

a surface tension 

Subscripts: 

c Jet core region 

g gas phase 

l,o inside and outside of a particular Jet region 

1 liquid phase 

IBTROOUCTZQN 

Liquid pumping can be obtained by utilizing 
the buoyant force of gas bubbles rising through the 
liquid as shown In fig. 1. This is a free convection 
type of process using two uaaiscible fluids that have 
a large difference in density (1) (2J . Because of tne 
large density difference made possible by using a gas 
and n liquid, large amounts of liquid pumping can be 
obtained. The large local density difference is in 
contrast to ordinary thermally driven single phase 
free convection where the density variations are usu- 
ally small. 
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Aa tat* resting application or bubble pumping is 

for ice prevention la lakes (J) f|) (£). Since the 
maximum density of water is at *° C, the rater in the 
bottom regioas of a quiescent late earn he as much as 
* c wanner then the surface regions when the air tem- 
perature i» below freezing. Ice fonmttion on the late 
tea be diminished or prevented by pumping the warmer 
N «tiots water toward the surface by teens of rising air 
babbles. Babble puraping is also used in fluidised 


nates, sad the rising ga s bub hie s . 
AXALYSIS 


* having m opening about 1-1/2 aa 
in diameter was placed at the bottom at a tank 0,6 m 
deep and 0.31 at stpare in horizontal cross section. 
Notion pictures were ted* at several rates at air flow 












From ref. fiiy. for bubble Reynolds maskers greeter 
than about 200, which would be reasonable for the large 
bubblas ia the present application, the £ is essential- 
ly constant and equals about 1,5. then in the core re- 
gion of local radius a c (x) in fig, 5, on the average 
l/{ K ♦ 1) of the vertical height is occupied by gas 
bubbles and U/(K * 1) is occupied by wake regions, 
this assumes for aatbanatical simplicity that the bub- 
bles are rising in a chain-bubble Idahlon ia which the 
bubbles are vertically spaced by the liquid wakes. 
Strictly speaking such a configuration exists only over 
a entail, range of flow rates. However, the results 
fro® this nodal say apply to other bubble regions. 

This is because, as discussed later, the turbulent ea- 
trainaaet of liquid by gas is minor compared with the 
«Btraiaaast by liquid. This tends to diminish the im- 
portance of the exact configuration of the rising gas. 
Since the gas weight flow rate is constant. 
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liquid being carried upward by the drag on the liq- 
uid arising free the buoyancy force of the bubbles. 
This model permits using different relations for the 
entrainment by the gas and by ike liquid regions. 

The flow field is assumed to be steady. Iso- 
thermal and fully turbulent. The liquid density is 
assumed constant, but the gas density varies with the 
local pressure in accordance with the perfect gas 
law; thus the bubbles expand as they rise through the 
liquid. The bubbles axe assumed to be sufficiently 
large that their drag is fully turbulent and hence 
they rise at a constant terminal velocity relative 
to the liquid. The local bubble velocity is as- 
sumed equal to the local liquid velocity plus the 
bubble terminal velocity, ft is assumed that the gas 
leaves the orifice with negligible upward momentum. 
After the gas is released, there is an adjustment re- 
gion within which the bubbles achieve their terminal 
velocity. For an installation in a lake cur river 
this region would be small compared with the total 
rise height and hence this region ia not taken into 
account. 

As shown in fig. 1 the rising jet How will turn 
at the liquid surface and move radially outward. At 
the surface ike vertical velocity component has to go 
to zero. however, to, incorporate this condition re- 
quires coupling the solutions for the jet and the 
turning zone which is a difficult analysis. Tim turn- 
ing zone is not accounted for here; it Is assumed that 
the jet continues to the surface as shown by the dot- 
ted extrapolation lines in fig. 1. Because of their 
strong vertical buoyancy force, the bubbles do not 
turn very much with the flow; hence they exert their 
pumping effect all the way to the surface. Also 
there is zero shear at the liquid surface which facil- 
itates the turning and tends to minimize the influence 
of the turning zone in the rising plume. These same 
types of assumptions are discus sec in ref, (igj for a 
free convection plume above a heated cylinder. Con- 
sequently an analysis without the effect of the turn- 
ing zone should yield reasonable vertical pumping 

See ..CastlMily 

For a single rising bubble, let the ratio of the 
wake volume to the bubble volume be a quantity I, 


Perfect Sea haw 

The pressure at height x above the nozzle is 
0,(1, - x). Then from the perfect gas law the gas den- 
sity at x is 
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where 

face. 


p is atjsospheric pressure at the liquid sur- 


ilottld Continuity 

The liquid continuity equation accounts for the 
liquid carried into the jet by turbulent entrainment. 
The entrained liquid goes into the liquid region sur- 
rounding the core, or into the bubble wakes which are 
growing as the rising bubbles expand. 
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The term on the right will be expressed later in terms 
of a turbulent entrainment relation. 

Gas and liquid Namentui 

In most applications involving air driven jets, 
the surrounding body of water is large; hence for 
simplicity in the present analysis, the surrounding 
liquid region will be assumed infinite. The analysis 
of a jet in a small container would be much more com- 
plicated because of liquid recirculation and the in- 
teraction with the container walls. 

The upward buoyancy force of the bubbles pro- 
duces a change in momentum (usually quite small) of 
the gas bubbles, a momentum change of the liquid in 
the bubble wakes, and a aoreentum change of the liquid 
in the region surrounding the bubbly core. This 
yields the momentum equation as, 
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"Top Hat" Distributions 

To Integrate eqs. (1), (3), and (4), "top hat" 
velocity profiles are assueed which have yielded good 
results for free convection pliaes (£), 


u g (x,r) - u (x) 0 < r < a c (x) 


u^x.r) 


{ Uj(x) a c (x)<r<a 
0 a(x) < r 


<»(*) 


(5a) 

(5b) 


(*) 

The u and u, are related by 
8 l 


u g (x) - Uj(x) ♦ u_ 


( 6 ) 


where u^ Is the terminal velocity of a single bubble 
rising in a large region of quiescent liquid. 

Insert eq. (5) into eqs. (1), (3), and (4) to 

obtain 
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Entrainment Function 

An expression is now needed for the entrain- 
ment on the right side of eq. (8). As discussed by 
Morton ( 12) for free convection plumes above fires, 
for a buoyant single phase jet the entrainment de- 
pends on an additional variable which is the density 
in the plume relative to the density of the surround- 
ing fluid. The entrainment is reduced when the plume 
density is smaller than the density of the surrounding 
fluid. The present situation involves two distinct 
phases that each maintain their separate identities; 
hence there will be two types of entrainment, liquid- 
liquid and gas-liquid. For the latter there is a lack 
of entrainment information when the density ratio is 
as small as that for air to water. As soon as same 
liquid is in motion, much of the entrainment Is by 
liquid entraining additional liquid which is reason- 
ably well understood. 

It was deduced in ref. (12) by using the infor- 
mation in ref. (13) , that for a single phase jet the 
usual jet entrainment coefficient E used when the 
jet and outer fluid have the same deftsity, should be 
modified when the jet and outer fluid densities are 
different. Ttye E Q is multiplied by the density ra- 
tio (pj/p j 1 ' 2 where p^ is the density in the jet 
and p iS outside the jet. For the present situa- 
tion the liqui .-liquid entrainment terms will there- 
fore have a unity ratio factor, while the gas -liquid 
terns will contain (p^/o,) 1 / 2 . Since the entrain- 
ment is so different for the two phases, the two- 
phase nature of the jet will be retained rather than 
trying to assign an average density to the entire jet. 

The entrainment also depends on the velocity of 
the Jet relative to its surroundings, and on the in- 


tarfaelal area between the Jet end the aurroundlng re- 
gion. These velocities and areas ere different for 
the gu and liquid portions. The liquid region moving 
at velocity u, as shown in fig. 3, entrains liquid 
from the quiescent fluid around it. The higher veloc- 
ity bubbly region Is moving at velocity u^ relative 
to Its aurroundlng liquid, and hence should enhance 
tha entrainment process. A well defined interfacial 
area bounding the bubbly region can only be obtained 
by Utilizing a simplified moael such as in fig. 3, and 
tblr area la used to obtain the magnitude of the en- 
trainment. However, if this entrainment were assumed 
to be retained In the core, the core would then be- 
came a liquid jet containing a bubbly core and the 
difficulty of defining the entrainment for this two- 
phase inner Jet would become the same as that for the 
original two-phase problem. Although the cemcept of a 
bubbly core Is used to provide a well defined inter- 
facial area, it is realized that the bubbly motion is 
actually more random. With these considerations in 
mind and in the absence of any better information on 
such a two -stage entrainment process, it is assumed 
that the vigorous action of the bubbly region in- 
creases the turbulence in the liquid surrounding the 
bubbles and thereby enhances the total entrainment In- 
to the moving region. This is in accord with the re- 
sults in ref. (13) that the entrainment is a function 
of the excess momentum flux in the jet. Thus in the 
present model, as the jet grows with increasing height 
above the source of gas the core will retain its iden- 
tity as being composed of only bubbles and ti. _>ir wakes. 
Hence the total entrainment will be taken as the en- 
trainment by the moving liquid outside the bubbly core, 
augmented by the entrainment effects of the liquid 
wakes and gas within ..us core, and is given by 
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solution for and a 

To obtain the amount of liquid being pumped by 
the rising gas, the liquid velocity and jet radius 
must be obtained. The liquid mass flow rate in the 
wake regions is known from the specified gas flow rate 
as KM_( o,/p ). The macs flow rate of liquid is then 

ft 8 

M l ■ ■ a c)v t + ^ (u) 

Since M in eq. (7) is a constant, eqs. (8) 
and (9) can be simplified by using eq. (7) to elim- 
inate some of the u Eq. (10) is then substitut- 
ed for the right siue of eq. (8). The result is 
the continuity and momentum equations in the form. 
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(13) 


The Ug can be eliminated In xerms of u^ by 
using eq. (6), and the p w Is given In tens of x 
by eq. (2). Fro* eq. (7) the a c can then be elim- 
inated in terms of u_ and x. Thus eqs. (12) and 
(13) are reduced to equations for the unknown Uj(x) 
and a(x). To obtain numerical results the differen- 
tiations were carried out analytically and then the 
equations combined to eliminate dUj/'dx or da/dx. 
The result was the following set of simultaneous dif- 
ferential equations that were solved by the Runge- 
Eutta method (14). 
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(15) 


At x ■ 0, a m a , and u ■ 0 so that u « u . Then 
c l 8 « 

the initial conditions at x - 0 - to begin the inte- 
gration are from eq. (7), 


RE8ULT8 

To compute the Induced liquid flaw, a nuaber of 
quantities must be specified. The results that follow 
are with regard to the prevention of ice in lakes by 
raising warm bottom water to the surface. For these 
conditions the bulk water is a few degrees K above 
fleeting so in what follows T - 275 K (495° R). The 
terminal velocity of a single bubble in undis- 
turbed fluid was obtained from the relation (15) (16) 


lbe(o, - cJT 1 ^ 4 

U.-1.53 - ' 

(16) 

l P l J 



which yielded ■ 0.252 m/sec (0.825 ft/sec) for air 
bubbles In water. 

A significant source of uncertainty is in 
specifying the entrainment coefficient E . In ref. 
(17) which Is concerned with the penetration of a 
c on d en s in g vapor jet into a liquid, a "top hat" ve- 
locity profile is used and a range of E. from 0.06 
to 0.12 la given as being found in the literature. 

For entrainment of various gas jets in still air with 
density ratios in the range d±/o 0 ■ 0.66 to 14.5, the 
work In ref. (13) yielded E q « o.O0 based on con- 
siderations of the excess momentm in the jet. The 
entrainment is proportional to a characteristic ve- 
locity in the jet. If a Gaussian velocity profile is 
used, the centerline velocity is used as the charac- 
teristic velocity which is larger than the average 
velocity used In the "top hat" profile. As a result 
In ref. (18) E « 0.082 is recommended for use with 
a Gaussian profile and E 0 » 0.116 for a "top hat" 
profile. In view of all these considerations it was 
decided to obtain two sets of calculations using E 0 => 
0.08 and 0.116. The results will he discussed in 
the next section. 

Equations (14) and (15) were integ? ed to 
x ■ L, the surface of the water, to obtain u, (L) 
and a(L). By using eqs. ( 6 ) and (7) a (L) Is found 
d the upward flow of liquid then obtained from 
eq. (11). In reality the upward flow will begin to 
turn as it approaches the surface, but the influence 
of the water surface was not accounted for here, as 
previously discussed. Figure 4 gives the mass flow 
rate of liquid being pumped to the surface as a func- 
tion of the mass flow rate of air, for the air being 
Introduced at various depths below the surface. In 
tens of the same quantities fig, 5 shows the radius 
of the jet region at x ■ L and fig. 6 gives the 
liquid velocity at x - L. Some of the results in 
fig. 5 have been cross plotted in fig. 7 to show the 
trend of the jet radius with orifice depth for fixed 
gee mass flow rates. 
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and u, (0) m 0. The latter condition however causes 
starting difficulty In the integration since it is in 
the denominator of a few terms, so a small value was 
used, Uj(0) ■ 0.001. The calculations were tested 
using smaller u (o) and no significant changes were 
found. 1 


DISCUSSION 


As shown by fig. 4 the liquid pumped to the 
surface Increases with the depth of the orifice below 
the surface and with the gas flow rate. For a fixed 
weight flow rate, the volume of the gas Introduced at 
the orifice decreases as the orifice depth increases. 
Hence it seems better to present the results in terms 
of gas weight flow rate than in terms of volume flow. 
For each orifice depth the results in fig, 4 lie fair- 
ly well along a straight line on a log -log plot. The 
slope decreases somewhat as L is Increased but the 
results vary essentially as M, « M’ . A cross plot 
of the results shows that for a given M the M 
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variation with L is also as a power function. 
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entralnmsnt tan omitted. As might be expected be- 
cause of the seen value of (o i^/p^) 1 ” the gee en- 
trainment term had less then a 1$ effect on the 
pwped liquid mass. The liquid wake ten had e sig- 
nificant effect for small x where the total liquid 
flow is small, as it helped initiate the entrainment 
process. The effect of this ten decreased as the 
total entrained flow became more substantial. For a 
depth of 1.524 m (5 ft) this ten Increased the flow 
in the range of 10 to 7?. as the gas flow ranged ti ram 
the smallest to the largest values on fig. 4. For a 
depth of 9.144 m (30 ft) the wake ten contributed 
about 5% to the pumped flow for any gas flow rate. 

The liquid wake ten on the left aide of eq. (9) 
is an especially important feature of the present an- 
alysis. For small x it provides inertia within the 
flow when the total entrained liquid is still small 
and prevents the buoyancy force from producing unre- 
alistically high liquid velocities near the origin of 
the numerical calculations. 

The bubble pumping process provides » fluctua- 
ting flow so that experimental measurements are dif- 
ficult. Seem comparisons will now be mads with the 
small amount of data available in the literature. 
These results will serve to emphasize that although 
there Is still appreciable uncertainty in the amount 
of flow being pumped, the trends with depth and gas 
flow rate have been approximately established. The 
most striking feature Is that the ratio of pwped 
liquid volume to gas volume is quits large. 

In fig. 13 of ref. (6j Kobus gives the .ratio of 
pumped liquid volume to the gas flow rate for various 
gas flow rates. There is a set of data for the local 
flow rate at a location 3.3 m above an orifice 4.5 m 
deep. These local flow rates were also obtained from 
the present computer calculations and a comparison is 
made in table I. For E Q » 0.8, the present theory 
predicts a flow that Is somewhat low at the high gat 
flow rates and somewhat high at the low flow rates. 

For E Q > 0.116 the agreement is good at the high 
air flow rates. Before further discussion the re- 
sults of ref. (£) should be included. 

In ref. (4) the following correlation is given 
from available data in terms of volume flow rates of 
air and water and depth of submergence ( in ref. (5) 
the exponent on is 1/2), 

Q m Const. 

Also in refs. (4) and (5) a discharge of 0.000472 m 3 / 
sec (0.0167 ft ’/sec) of air at a submergence of _ 
1.68 m (5.5 ft) is reported to raise about 0.0567 m V 
sec (2 ft 3 /sec) of water to the surface. The present 
calculation yields a pumping of about 0.0373 mY'iec 
(1.33 ft 3 / sec) of water when E 0 « 0.08 is used, and 
0.0585 m 3 /sec (2.07 ftVsec) when E - 0.116. The 
experimental values were inserted into eq. (18) to 
evaluate the constant and the formula then used to 
calculate local values for the conditions In table I. 
The flow rates are seen to be about twice those given 
in ref. (£) . 

It is noted from table ! that the pumped vol- 
ume of water per unit volume of air decreases as the 
air flow rate increases, so that a more efficient 
pumping system is obtained at low flow rates. In 
ref. (£) it was found that the volume pumping rate 
decreased as the air flow to the -0.4 power, in 
ref. (Jj it is the -0.33 power, in ref. (5) it is the 
). 5 power and the present eq. (17) gives approximat- 
ely the -0.6 power. It is felt that additional sets 
of data are needed before refinements can be made in 


the theory. 

To provide a little more Information on bubble 
pi eying consider briefly some results reported in ref. 
(£) from an approximate analysis made in 1946 by 
Konovalov of the pumping in a geometry that is two 
dimensional in rectangular coordinates as produced by 
air released from a submerged perforated pipe. The 
empirical constants in the theory were evaluated from 
laboratory observations for depths of the perforated 
pipe up to about 1 a and air flow rate; up to 0.00138 
mveec par meter of pipe length. There were also some 
limited full scale tests run at depths up to 10 m. 

The results for orifice depths greeter than 1 m are 
given by the correlation 

^ - 0.75 £(10 4 L)L 2 In (l + ( 19 ) 

where and are in m 3 /sec and 1 is in meters. 

Because of the difference in geometry precise compar- 
isons cannot be made, but it is interesting to see how 
the predictions of liquid volume pumped per unit gas 
volume Qj/Q Caspars for the release from an orifice 
as obtained Tram eq. (17), and from a meter length of 
perforated pipe. eq. (19). Results are given in ta- 
ble II and it is seen that the pumping rates are of 
the same magnitude. The was assumed to be at the 
discharge location, although this was not clearly spec- 
ified in the reference. 

In ref. Qg) it is mentioned that the bubbly re- 
gion of the axi symmetric Jet is approximately con- 
tained within a cone having a total included angle of 
12°. Although the present mathematical model consid- 
ers ell the gas to be in the form of large bubbles 
contained in the core of the Jet. In the physical case 
there ere s mall bubbles that break off from the large 
ones as shown in fig. 2. The turbulent motion diffuse* 
the small bubbles within the flow so that they prob- 
ably extend throughout most of the entrained region. 
Hence the totel bubbly region should give an indica- 
ted the extent of the jet region. Using a total 
inc.ulei cone angle of 12°. the jet radius at the sur- 
face is 

a (x • L) ■ L tan 6° ■ 0.105 L (20) 

The photographs in fig. 2 indicate an increase of cone 
angle with gas flow rate. In fig. 2(a) the cone angle 
is about 14°, while in fig. 2(d) it is about 19°. 
Equation (20) is plotted in fig. 7 and compares reas- 
onably well with calculated results for an entrainment 
coefficient of E- « 0.08. For a larger total includ- 
ed angle of say 18°, which is characteristic of figs. 
2(c) and (d), eq. (20) becomes a(L) » 0.158 L. As 
shown by fig. 7, this provides reasonable agreement 
with the values computed with E Q - 0.116. 

CONCLUSIONS 

A tractable mathematical model was formulated for 
computing the liquid carried upward in an axisymmetric 
jet driven by a rising stream of gas bubbles. The 
model geometry was based on observations of air rising 
through water 0.6 m deep. Liquid is carried into the 
Jet by turbulent entrainment and a difficulty in the 
•analysis is in specifying the proper value of the en- 
trainment coefficient. Calculations were made for two 
values of the entrainment coefficient within the range 
given in the literature. The large buoyancy resulting 
from the large density difference between the gas and 
liquid produced considerable liquid movement compared 
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with ordinary liquid free convection induced by 
tbemel means, where density difference* are usually 
on the order of a few percent. The present result* 
indicated that for air rising through water the rate 
of liquid transport to the water surface varies as 
the 0.4 power of the air mass flow rate and the 1.4 
power of the submergence depth of the air orifice. 

More experimental data is needed before the theory 
can be further refined. 
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TABLE I. - VOLUMETRIC PUMPING RATES IN AN aXI- 
SYMMETRIC JET AT A LOCATION 3.3 M ABOVE ORIFICE SUB- 


MERGED 4.5 M BELOW SURFACE 


Air flow 
at standard 
conditions. 
Q g . * Veec 

Volume water manned Ql 
Volume air Qg 

Kobus (6) 

Baines (±) 
( extrap- 
olated) 

Present c< 

e 0 - o.oe 

ilculations 
E 0 . 0.116 

0.00040 

175 

349 

217 

327 

.00130 

100 

236 

108 

163 

.00255 

86 

188 

73 

no 

.00420 

67 

159 

54 

81.5 

.00620 

62 

140 

42 

63 


TABLE II. - VOLUMETRIC FUMPING RATES FOR SUBMERGED 


AXI SYMMETRIC JET OR PERFORATED PI IE (Q g IS VOLUME 


FLOW AT ORIGIN OF GAS DISCHARGE) 


Pipe 

depth, 

L 

m 

*8 

m 3 /sec 

*1 

Axisymme 
from e 

E q - 0.08 

/Q g 

trie jet. 
q. (17) 

E q . 0.116 

“A 

Perforated pipe 
1 m long, 
eq. (19) 

2 

0.0001 

284 

440 

717 


.001 

71.4 

111 

154 


.01 

17.9 

27.8 

33.3 

3 

0.0001 

540 

836 

1090 


.001 

136 

211 

234 


.01 

34.1 

52.8 

50.6 

6 

0,0001 

1490 

2310 

2250 


.001 

374 

580 

485 


.01 

94.0 

146 

105 

12 

0.0001 

4460 

6910 

4720 


.001 

1120 

1740 

1020 


.01 

281 

435 

219 



